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In standard general-relativistic cosmology, for fluids with a linear equation of state (EoS) P = wp 
or scalar fields, the high isotropy of the universe requires special initial conditions: singularities are 
velocity dominated and anisotropic in general. In brane world effective 4-dimensional cosmological 
models an effective term, quadratic in the energy density, appears in the evolution equations, which 
has been shown to be responsible for the suppression of anisotropy and inhomogeneities at the 
singularity under reasonable assumptions. Thus in the brane world isotropy is generically built 
in, and singularities are matter dominated. There is no reason why the effective EoS of matter 
should be linear at the highest energies, and an effective non-linear EoS may describe dark energy 
or unified dark matter (Paper 10). In view of this, here we investigate the effects of a quadratic 
EoS in homogenous and inhomogeneous anisotropic cosmological models in general relativity, in 
order to understand if in this standard context the quadratic EoS can isotropize the universe at 
early times. With respect to Paper iQ, here we use the simplified EoS P — ctp + p 2 /p c , which 
still allows for an effective cosmological constant and phantom behavior, and is general enough to 
analyze the dynamics at high energies. We first study homogenous and anisotropic Bianchi I and V 
models, focusing on singularities. Using dynamical systems methods, we find the fixed points of the 
system and study their stability. We find that models with standard non-phantom behavior are in 
general asymptotic in the past to an isotropic fixed point IS, i.e. in these models even an arbitrarily 
large anisotropy is suppressed in the past: the singularity is matter dominated. Using covariant and 
gauge invariant variables, we then study linear anisotropic and inhomogeneous perturbations about 
the homogenous and isotropic spatially flat models with a quadratic EoS. We find that, in the large 
scale limit, all perturbations decay asymptotically in the past, indicating that the isotropic fixed 
point IS is the general asymptotic past attractor for non phantom inhomogeneous models with a 
quadratic EoS. 

PACS numbers: 98.80.-k, 95.35.+d, 95.36.-fx, 98.80.Jk 



I. INTRODUCTION 

Observations of large scale structure (LSS) and cos- 
mic microwave background radiation (CMBR) 0, H Q 
indicate that the universe is highly homogenous and 
isotropic on large scales. Under standard assumptions 
on the matter content, in the standard general relativistic 
cosmological models isotropy is a special feature, requir- 
ing a high degree of fine tuning (a special set of initial 
conditions) in order to reproduce the observed universe. 
Specifically, in the set of spatially homogeneous cosmo- 
logical models with a fluid which satisfy the energy condi- 
tions p > \P\ and P > 0, those which approach isotropy 
at late times is of measure zero . In these models the 
anisotropy dominates the dynamics in the neighborhood 
of the singularity, while the effect of matter is negligible: 
thus the singularity can be said to be velocity dominated 
0- In general anisotropic models do not isotropize suf- 
ficiently as they evolve into the future (models do not 
evolve towards a Friedmann-Robertson- Walker universe 
rapidly enough). This problem is known as the isotropy 
problem and can be solved by inflation. Many investi- 
gations on this issue show that a large class of models 
evolve towards a homogenous and isotropic model when 



under the influence of inflationary matter. This means 
that even models with initially high levels of anisotropy 
can evolve into the universe we observe today. A large 
body of work exists relating specifically to the study of 
spatially homogenous and anisotropic Bianchi cosmolog- 
ical models. In a classical paper, Wald showed that ini- 
tially expanding Bianchi models with a cosmological con- 
stant asymptotically approach a de Sitter state (hence a 
spatially homogenous and isotropic model), except for 
the Bianchi type IX mixmaster models which may re- 
collapse 8J. In the case of inhomogeneous models it is 
not so clear cut and in order for inflation to begin, suf- 
ficiently homogenous initial conditions are required (see 
e.g. 0,0|). As far as singularities in general relativity 
are concerned, it is widely accepted Hlj that the Belin- 
skii, Khalatnikov and Lifshitz (BKL)[l2L Il3| picture is 
correct: for fluids with linear EoS P = wp (w < 1) or 
scalar fields 1 , when the singularity of generic inhomoge- 
neous anisotropic cosmological models is approached, the 



Massless scalar fields are equivalent to a stiff fluid P = p (w = 1) 
and are special in this respect, see 1 1 -J and references therein. 
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dynamics is asymptotically that of the mixmaster model, 
the singularity is thus velocity dominated and the effect 
of matter is negligible 0] . 

More recently various theories of gravity that go be- 
yond general relativity have been considered. For in- 
stance, high energy modifications of gravity could come 
from extra dimensions (as re quir ed in string theory). In 
the brane world 

Ella III Hi 

scenario the extra di- 
mensions produce a term quadratic in the energy density 
in the 4-dimensional effective energy-momentum tensor, 
i.e. this quadratic term appears in the 4-d evolution and 
constraint equations on the brane. Under the reasonable 
assumption of neglecting 5-d Weyl tensor contributions 
in these equations, this quadratic term has the very in- 
teresting effect of suppressing anisotropy at early enough 
times, as the singularity is approached. In the case of a 
Bianchi I brane- world cosmology containing a scalar field 
with a large kinetic term the initial expansion is quasi- 
isotropic |l9j . In this scenario, the initial anisotropy pro- 
vides increased damping in the scalar field equation of 
motion, thus resulting in an extended period of infla- 
tion. The rapid decay of the anisotropy results in the 
class of initial conditions allowed by observations to be 
greatly increased. Under the same assumptions, Bianchi 
I and Bianchi V brane-world cosmological models con- 
taining standard cosmological fluids with linear equation 
of state (EoS) (P/p = w > 0) also behave in a similar 
manner |20j. The same is also true for more general ho- 
mogeneous mod els 2lll22| and even some inhomogeneous 
exact solutions 23] . Finally, within the limitations of 
a perturbative treatment, the quadratic-term-dominated 
isotropic brane-world models have been shown to be local 
past attractors in the lar ger s tate space of inhomogeneous 
and anisotropic models |24l [25| . More precisely, again 
assuming that the 5-d Weyl tensor contribution to the 
brane can be neglected, large-scale perturbations of the 
isotropic models vanish in the past, approaching the sin- 
gularity. Large scale perturbations A ^S> are the only 
type relevant to this kind of analysis for non-inflationary 
matter, because for any given A there is always an early 
enough time such that A 3> H^ 1 is satisfied 2 . Thus in the 
brane scenario the observed high isotropy of the universe 
is built in, i.e. the natural outcome of generic initial con- 
ditions, unlike in general relativity where general cosmo- 
logical models with a linear EoS are highly anisotropic in 
the past, approaching the singularity a-la BKL, i.e. with 
mixmaster dynamics ^2 13 13 • 



2 We are referring here to the well known fact that during a non- 
inflationary (i.e. decelerated) phase a < a given scale much 
larger than the Hubble horizon A 2> H~ 1 at a given time will 
"enter the horizon" ( A = H -1 ) at a later time. In practice, 
neglecting small scales perturbations is equivalent to neglect gra- 
dients and thus the Laplace operator term in the evolution equa- 
tions, which then become ordinary differential equations. In the 
perturb ative context, this is analogous to the BKL approxima- 
tion mm 



It would therefore be interesting to see if the behavior 
of the anisotropy at the singularity found in the brane 
scenario is recreated in the general relativistic context 
if we consider a quadratic term in the EoS. That is, 
do we have a scenario where, in general relativity, the 
generic asymptotic attractor in the past is an isotropic 
model, thanks to a non-linear EoS. This is the question 
we aim to investigate, for the specific case of the EoS 
P = ap + p 2 1 p c . This is a simplified version of that used 
in Paper I [l| , where the focus was on using a quadratic 
EoS to represent a dark energy component or unified dark 
matter. It is general enough to analyze the high energy 
dynamics we are interested here, while still allowing for 
an effective cosmological constant and phantom behav- 
ior. 

The paper is organized as follows. In section II we out- 
line the setup of the models under investigation, defining 
three classes, one of which represents a fluid with phan- 
tom behavior. Using dynamical system methods [26ll2^| . 
in section III we consider the anisotropic Bianchi I and 
V class of homogenous models containing a perfect fluid 
with quadratic EoS, showing how models with standard 
non-phantom matter are asymptotic in the past to an 
isotropic fixed point IS. Using covariant and gauge in- 
variant variables [2^, |2^, |3(j , in section IV we give the 
linearized evolution and constraint equations for generic 
inhomogeneous and anisotropic perturbations about a 
flat Friedmann model with a quadratic EoS. We solve the 
equations in the large scale limit A 3> H _1 and show that 
the perturbations vanish as the singularity is approached, 
indicating that the isotropic model IS is the generic past 
attractor. We then finish with some concluding remarks 
and discussion in section V. 



II. THE MODEL 

We begin with a summary of the most relevant points 
that can be derived mostly from an analysis of energy 
conservation alone. More details can be found in Section 
II in Paper I [l| . 

The energy conservation equation for a cosmological 
model containing a perfect fluid is: 

p = -3H(p + P), (1) 

where p is the energy density, P the isotropic pressure 
and H is the Hubble expansion scalar. An average scale 
factor a is defined, at least locally, by a/a = H. Through 
this relation solutions of Eq. (JTJ can be found in terms 
of a, but first an EoS must be specified. Usually this is 
taken to be a linear function of the energy density of the 
form P = wp, with w = constant. In this case Eq. ifTJl 
gives: 



In the case of the homogenous and anisotropic Bianchi 
I cosmological models, the amplitude of the shear scales 
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(3) 



This example shows how the shear can decay rapidly in 
an expanding model, growing larger for smaller a. If the 
dynamics of the model is such that a — ► 0, both p and 
a blow up, and the model has a physical singularity. In 
the case w < —1, the so called "phantom" fluid |3jJ, the 
energy density counter intuitively decays (grows) in the 
past (future) and therefore any shear component always 
comes to dominate at the initial "Big Bang" singularity. 
In the case when — 1 < w < 1 the energy density decays 
during expansion, but at a slower rate than the shear, 
so that once again the shear component dominates at 
the singularity. In approaching these singularities matter 
with linear EoS w < 1 is un-influential, therefore these 
singularities may be termed velocity dominated 0. It 
is only in the case when w > 1 ("super-stiff' fluid) that 
the energy density is dominant in the past and therefore 
the initial singularity is isotropic. This indicates that 
as the initial singularity is approached in the standard 
homogenous model with w < 1, any residual anisotropy 
will generally come to dominate near the singularity. The 
anisotropy is a generic feature of these models and the 
initial singularity is in general anisotropic and velocity 
dominated. 

As discussed earlier, motivated by the results in the 
brane-world models, we aim to investigate the effects 
of introducing a non-linear EoS on the behavior of 
anisotropic models. The EoS we propose is quadratic in 
the energy density (we neglect the constant pressure term 
introduced in Paper 10 as in general this only modifies 
late time behavior): 



P = ap + — . 



(4) 



The discrete parameter e denotes the sign of the 
quadratic term, e £ { — 1,1}, and p c > sets the en- 
ergy scale above which this term becomes relevant. The 
energy conservation equation can now be integrated 
to give the scaling of the energy density as a function of 
the scale factor (except for a = —1). 
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A(a + l)p c 

ffl 3(a+l) _ ej 4' 
3(a+l) 

(a + l)p c + ep ' 



(5) 
(6) 



where p Ql a a represent the energy density and scale factor 
at an arbitrary time t . This is valid for all values of e, 
p c and a (a ^ —1). Defining 



PA = -e(l + a)p c 



(7) 



the fluid admits an effective positive cosmological con- 
stant point, pa > 0, only if e(l + a) < 0. As shown 



in Paper I [| , in the case of e = — 1 the fluid generally 
behaves either in a phantom manner or is asymptotic to 
a finite energy density in the past, the effective cosmo- 
logical constant p\, for all open and flat models. In this 
scenario, the initial singularity is going to be dominated 
by anisotropic expansion as seen in the standard linear 
EoS case, due to the sub-dominance of the energy density 
at early times; therefore we do not consider the e = — 1 
case further. 

In the case of e = +1, the expression (JSJ for the energy 
density p(a) is conveniently rewritten in three different 
ways, defining a+ — | J 4| 1 / 3 ( Q + 1 ) and assuming p > 0. 



A: (1 + of) > 0, pa < 0, 



(1 + a)p c 



(i) 



3(l+ct) 



(8) 
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In this case a* < a < oo, with oo > p > 0. The past 
singularity is similar to certain future singularities found 
in dark energy models [32, H3 and is referred to as Type 
III singularities: 

• For t — > t*, a — > a*, p — > oo and |P| — > oo 

Where and a* are constants with a* ^ (we reach a» 
at a finite £*). The main difference is that in our case 
the singularity occurs in the past. At late times the fluid 
behaves in a similar manner to a fluid with a linear EoS 
and w = a. 



B: (1 + a) < 0, < p A < p, 



PA 



1 



(i) 



3(l+a) ' 



(9) 



In this case a* < a < oo, with oo > p > p A . As in 
case A, the fluid approaches a Type III singularity in 
the past, but now approaches an effective cosmological 
constant at late times. 



C: (1 + a) < 0, < p < p Al 



P = 



Pa 



1+ ^ 



3(l+a) 



(10) 



In this case < a < oo, with < p < p\. The fluid 
behaves in a phantom manner but tends to an effective 
cosmological constant at late times. 

In the cases A and B, as we approach the past singu- 
larity the energy density approaches infinity, in fact the 
singularity occurs at a finite scale factor, that is for the 
energy density one has p — > oo as a — * a*. The cosmolog- 
ical model emerges at a finite scale factor a*, with infi- 
nite energy density and infinite isotropic pressure P. The 
model emerges from a curvature singularity, infinitely de- 
celerates and asymptotically approaches a standard flat 
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Friedmann model with a linear EoS (case A) or a de 
Sitter model (case B). In the case C, the fluid behaves 
in a phantom manner and approaches p\ at late times. 
In this case the energy density is sub-dominant at early 
times and the anisotropy will dominate as we will show 
in the proceeding sections. In cases A and B, we have 
the possibility that the energy density could be the dom- 
inant component in the past, forcing the singularity to 
be isotropic. This would result in isotropic expansion at 
early times as seen in the brane-world model, with the 
fluid also approaching a non super-stiff EoS at late times. 
We will now proceed with the analysis of anisotropic cos- 
mological models containing a fluid with a quadratic EoS. 

III. DYNAMICS OF BIANCHI MODELS 

We now study the dynamics of homogeneous and 
anisotropic Bianchi cosmological models. In particular 
we investigate the Bianchi I and V classes of models that 
contain the flat and open Friedmann models. The matter 
content will be modeled by a non-tilted (that is the fluid 
flow remains orthogonal to the group orbits) perfect fluid 
with a quadratic EoS. 

We make use of the orthonormal frame formalism 
which results in a system of first-order evolution equa- 
tions. We have adopted a notation consistent with that 
used in [2(| (Greek indices run from a = 1, ... ,3 and 
Latin indices run from i = 0, ... ,3). The orthonormal 
frame formalism is a 1+3 decomposition of the Einstein 
field equations into evolution and constraint equations. 
The decomposition is carried out relative to the time- 
like vector field eo of a given orthonormal frame {e^}. In 
the case of models which admit an isometry group, eo 
is chosen to be the normal (u) to the space-like group 
orbits. The Bianchi models are such cosmological mod- 
els, they admit a three dimensional group of isometries 
which act simply and transitively on the space-like hyper- 
surfaces (surfaces of homogeneity). The orthonormal 
frame {u, e Q }, is such that: 

u u = —1, u • e Q — 0, e a • = 5 a p. (11) 

The dynamics can then be described in terms of the spa- 
tial commutation functions, the kinematical quantities 
and the matter variables. The spatial commutation func- 
tions, 7 Q /3/i, are defined by the commutation relations 
between the spatial basis vectors such that: 

[e a ,ep] = j^ape^. (12) 

The spatial commutation functions can be decomposed 
into a 2-index symmetric object and a 1-index object 
which determine the connections on the space-like 3- 
surfaces. 

l a ^ = e^ y n av + a p S ll a -a^5 p a . (13) 

The kinematical quantities include the Hubble scalar H , 
the shear tensor a a p and the angular velocity fl a , relative 



to a Fermi-propagated spatial frame, of the spatial frame 
{e Q } along the time-like vector u. The relevant matter 
variables are the energy density p and the isotropic pres- 
sure P, while the energy flux and the anisotropic pressure 
are set to zero as we are considering a perfect fluid. The 
decomposed Einstein field equations for a general Bianchi 
model containing a perfect fluid then take the form: 

H = -H 2 -\n 2 -\{p + 2,P), (14) 
3 6 

<j Q/ 3 = —3H<j a p — ^ S a/ 3 

+2n v (s^ a a 0fi + e"V aM ), (15) 

p = 3H 2 - a 2 + i (3) i?, (16) 

= 3a a % - e/'V/n^. (17) 

The trace- free spatial Ricci tensor ( 3 S a p = 3 R a p — 
| 3 R8 a p) and the spatial curvature scalar ( 3 R = 3 R^) 
can be given in terms of the spatial commutation func- 
tions: 

3 S a/ 3 = b a p — -(b^Sap 

^a^e^anp^ + e^pn^), (18) 

3 R = -iftp" - 6a p a". (19) 

The quantity b a p, can be expressed in terms of the spatial 
commutation functions and has the form: 

b a p = 2n c /n llfs - n^n a p. (20) 

The Jacobi identity for vector fields when applied to the 
spatial frame vectors e Q , in conjunction with the Einstein 
field equations give a set of evolution and constraints 
equations: 

n a p = -Hn a p + 2(er f /«-/3 AI + o^n apl ) 

+2Q,{e^ a ap tl +s^pa atl ), (21) 

a a = -Ha a - a a f3 ap + Ea^a^Vly, (22) 

= njap. (23) 

The evolution equation for the energy density results 
from the contracted Bianchi identities, or more directly 
from the conservation equations, and has the form QJ. 
The latter with the above system of equations reduce 
to the evolution and constraint equations for given non- 
tilted perfect fluid Bianchi model by making further re- 
strictions on the spatial commutation functions (n a p,a a ) 
and kinematical quantities (a a p,fl a ). The Bianchi I and 
V models are considered separately in the following sub- 
sections. 



A. Bianchi I cosmologies 

The Bianchi I models are a subclass of the Bianchi 
Class A models (a a = 0) . These models are homogeneous 
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and anisotropic cosmological models containing the flat 
Fricdmann (F) model. In the case of Bianchi I, the group 
invariant orthonormal frame {u,e Q } can be further spe- 
cialized such that: 



0, 
0. 



(24) 
(25) 



This in turn means that the three curvature of the spatial 
hyper-surfaces ifTTfl H$jl is zero: 



•a/3 



0. 



(26) 



The Einstein field equations l|15l I17fl and the Jacobi iden- 
tity lf2*T)l imply: 



a a j3 = diag (ai i, 022, 033), il a = 0. 



(27) 



The analysis of the system of equations can be simplified 
by introducing the following set of normalized dimension- 
less variables: 



Pe- 
er 



H 



fp~ c t 



(28) 



where a = Oapo^ jl. In addition, it is customary to 
introduce 



H y 



(29) 



as a dynamically dimensionless measure of anisotropy: 
we can refer to a cosmological model as isotropic in some 
limit if, in that limit, £ — ► 0. In particular, if there is 
a singularity for some a s (possibly zero), we refer to the 
singularity as isotropic if £ — > in the limit a — > a s . 
The above normalization to p c is useful because the en- 
ergy scale p c itself disappears from the dynamical system: 
in terms of the dimensionless variables above, the com- 
plete system of equations for Bianchi I models is explicitly 
made scale invariant, and is: 



V 
x' 

y' 

z 



x z 
3 + T' 

-3y ((a + l)x + x 2 ) , 
2 1 

-y 
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((3a + l)x + 3x 2 ) 



-3yz. 



(30) 
(31) 
(32) 
(33) 



The primes denote differentiation with respect to 77, the 
normalized time variable. The variable x is the normal- 
ized energy density, y is the normalized Hubble scalar 
and z is the normalized shear. The effective cosmological 
constant point is now given by x\ = — (1 + a). We only 
consider the region of the state space for which the en- 
ergy density remains positive [x > 0) and models which 



are expanding (y > 0). The evolution equation for the 
shear can be integrated to give: 



(34) 



The shear tends to grow in the past and rapidly decays 
at late times. Using Eqs. ©, ©, (JD) and C2} we 
can see how the shear scales with respect to the energy 
density (by substituting for the scale factor in each case) 
for each of the three sub-cases of the fluid. 



Case A: (1 + a) > 0, pa < 0, 



a + 1 +x 



(35) 



In this case, at early times the energy density diverges, 
x — > 00 and the shear goes to a constant value, z — > z a . 
The initial singularity is dominated by the energy 
density and is forced to be isotropic, as S — > 0. At late 
times, x — > and z — * 0. 



Case B: (1 + a) < 0, < p A < p, 



x — x\ 



(36) 



In this case, at early times the energy density diverges, 
x — » 00; the shear goes to a constant, z — > z a . The initial 
singularity is again dominated by the energy density 
and is forced to be isotropic, with S — > 0. At late times 
the energy density approaches a constant value, x — > a; a 
and z — + 0. 



Case C: (1 + a) < 0, < p < pa, 

x A - x \ T^+r 



(37) 



In the final case, the fluid behaves in a phantom manner. 
At early times x —> and z — » 00. The initial singularity 
is now dominated by the shear and is anisotropic. At 
late times the energy density approaches a constant 
value, x — ► x\ and z — * 0. 

In cases A and B we expect the generic past attrac- 
tor to be an isotropic singularity, while at late times we 
expect either a Minkowski or de Sitter model to be the 
attractor. In case C, the fluid behaves in a phantom 
manner and we expect models to be asymptotic to a 
anisotropic singularity in the past and a de Sitter model 
in the future. We now carry out a dynamical systems 
analysis of the above system. The system can be reduced 
from a 3D system to a 2D system by substituting for a 
variable using the generalized Friedmann equation l|3U|l . 
We choose to remove the normalized Hubble scalar y, as- 
suming expansion, y > 0. The resulting reduced planar 
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dynamical system is then: 



x 
z' 



) ((a + l)x + x 2 ) , (38) 
-y/3(x + z 2 )z. (39) 



The above system of equations is of the form = 
fi(uj). Since this system is autonomous, trajectories 
in state space connect the fixed/equilibrium points of 
the system (uj. ), which satisfy the system of equations 
fi{ u j,o) — 0. The stability nature of the fixed points 
can be found by carrying out a linear stability analysis, 
see Section II. B in Paper I [lj for a summary, and e.g. 
[2^,123 for more details. The fixed points of the reduced 
planar system JSSJl-(JS3J|, the corresponding eigenvalues 
and there existence conditions (the energy density re- 
mains positive (x > 0) and all three variables are real 
(x, y, z G M)) are given in Table |I] 



TABLE I: Location, eigenvalues and existence conditions (x > 

and x, y, z £ R) for the fixed points of the reduced Bianchi 

1 planar system 1381 - 1390 . To simplify we use r = — (a + 1). 



Name 


x, z 


Eigenvalues 


Existence 


M 


(0, 0) 


0, 


a£l 


dS 


(r, 0) 




a < — 1 



The first point represents the vacuum static Minkowski 
model. Eigenvalues of the linearized system are zero, 
hence the linearization theorem |27| doesn't apply at this 
point and we resort to numerical analysis to determine 
the stability character. In the case a > — 1, this is the 
only finite valued fixed point present. This is clearly 
seen in Fig. [Ja), corresponding to case A, where the 
black lines represent separatrices, the dots represent fixed 
points and the thin grey lines are example trajectories. 
The Minkowski point M is clearly the future attractor 
for all trajectories in the region. The behavior in the 
past is not clear from this diagram (but trajectories do 
appear to evolve to a constant z). The horizontal black 
line represents the Bianchi I Kasner vacuum models, with 
line element: 

ds 2 = ~dt 2 + t 2pi dx\ + t 2p2 dx\ + t 2p3 dx\ . (40) 
The parameters p a satisfy the following constraints: 



Although these models are asymptotic to Minkowski 
at late times, as seen in Fig. ^a), they are always 
anisotropic, in the sense that £ = v3 at all times. 
The Kasner singularity is a paradigmatic, characteristic 
of many models; the more general Bianchi IX singular- 
ity is approximated by a series of jumps from a Kasner 
phase to another, see e.g. 0,H(|. The vertical black line 
represents the expanding flat Friedmann model with a 
quadratic EoS. In the case of a spatially flat model as in 
this case, the Friedmann equation can be integrated to 
give the dependence of the scale factor on time: 



[3(a- 



arctan 



I)] 3 / 2 



^/(a/a*) 3(Q+1) -l 



(43) 



The scale factor starts at finite size a+ (we have fixed 
constants so that this corresponds to 77* = 0), expands 
and asymptotically approaches the behavior of a model 
with the standard linear EoS, 



a w a+(rj - T] ) 



2/3(c 



(44) 



For a < — 1, Fig. ^b), there is a second fixed point 
dS representing the generalized flat expanding de Sitter 
model. This point has attractor stability character, i.e. in 
this case the generalized de Sitter model is the generic fu- 
ture attractor. The new horizontal black line in Fig.^b) 
is a separatrix representing models with ia = — (a + 1); 
it is the dividing line between the phantom (x < Xa) 
and standard behavior (x > xa)- The trajectories be- 
low this line represent phantom models (corresponding 
to case C), that is the energy density increases as the 
model expands. The models in this region appear to 
be dominated by shear in the past (x — > 0, z — > 00 as 
rj — > —00) and asymptotically approach a de Sitter model 
in the future [x — > xa), z — > as rj — > 00). The scale 
factor for an expanding flat de Sitter model has the form: 



(45) 



The trajectories above the separatrix (corresponding to 
case B) also asymptotically approach the generalized de 
Sitter model in the future. 

The behavior in the past is not clear from the state 
space in Fig. ^ In order to analyze the dynamics at 
infinity, we introduce compactified Poincare variables X,; 
these can be expressed in terms of the standard variables 
Xj as: 



(41) 



71=1 



n=l 



and the average scale factor for expanding models of this 
type is given by: 



a = a (r/ - r/ a ) 1/3 . 



(42) 



X, 



The relationship can be inverted to give: 

Xi 



(46) 



(47) 



2 4 6 8 10 ' 2 ' ' 4 ' 6 ' 8 ' 10 

z z 

FIG. 1: The state space for the Bianchi I reduced planar system 1381 1- 139( 1. showing fixed points at finite values of the energy 
density and shear, (a) Left panel: the state space for a > — 1; the only fixed point present represents a Minkowski model, (b) 
Right panel: the state space for a < — I; in this case there is a new fixed point dS which represents a flat expanding de Sitter 
model with pa = — (1 + a)p c - The separatrix represents models with p = pA- Models below this line have phantom behavior. 



The new variables allow one to analyze the dynamics at 
infinity more easily, as Xi — > oo =>• Xi — > 1 . The evolution 
equations for the new variables are then given by: 



X> = 



\ 



E 

k=l 



X, 



k X k 



(48) 



where the x[ represent the evolution equations for the 
original variables. We now study the compactified and 
reduced system {X, Z); the fixed points of the system 
are given in Table ITU 

The first two fixed points (M and dS) represent the 
Minkowski and generalized de Sitter points as in the non- 
compactified case and retain the same stability character. 
The third point (IS) represents the isotropic singularity, 
which is matter dominated {X — > 1, or x — ► oo). This 
fixed point always has repeller stability and is the generic 
past attractor. The final fixed point (KS) represents the 
anisotropic Kasner singularity (Z — > 1, or z — > oo). This 
fixed point has saddle stability. The full state space for 
the system can now be reconstructed, and is depicted in 
Fig. |21 where as before the black lines represent separa- 
trices, the dots represent fixed points and the thin grey 
lines are example trajectories. 



In the case a > — 1, Fig. |2a) (this corresponds to case 
A, i.e. this is the compactified version of Fig. ^a)), 




TABLE II: Location and existence conditions (X > and 
X, Zel) for the compactified planar Bianchi I system. To 
simplify we use t = — (a + 1). The * is shown when the point 
does not exist. 

Name X, Z Existence Stability (a < — 1, a > — 1) 
(0, 0) a G R Saddle, Attractor 

() j a < — I Attractor, 

a£l Repeller, Repeller 

a G R Saddle, Saddle 



there are three fixed points, M, IS and KS. The hori- 
zontal black line X = represents the anisotropic Kas- 
ner model. The vertical black line Z — represents the 
flat Friedmann model containing a fluid with a quadratic 
EoS. The Minkowski point M is the future attractor for 
all trajectories in the region. The isotropic singularity IS 
is the generic past attractor and the anisotropic Kasner 
singularity KS has saddle stability character. The set of 
initial conditions for which anisotropy dominates in the 
past are special, corresponding to vacuum models X = 0, 
as opposed to being generic as in the case of the linear 
EoS. 
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FIG. 2: The compactified state space (X, Z) for the Bianchi I reduced system, including fixed points at infinity. These 
represent the isotropic singularity (IS) and the Kasner anisotropic singularity (KS). (a) Left panel: the state space for a > — 1, 
corresponding to case A, Fig. ITTa - ), (b) Right panel: the state space for a < — 1. The separatrix corresponds to that of Fig. 
1(b), with p = pa- Above this line we have models of type B. Models below this line are phantom, of type C, and Kasner-like 
at the singularity KS. Clearly the non-phantom models of type A and B are asymptotic to the IS fixed point in the past. 



In the case a < — 1, Fig. |2Jb), there are four fixed 
points, M, dS, IS and KS. The state space is split into 
two regions by the separatrix joining the points dS and 
KS and given by X = X\, where: 



X A = T 



1-Z 2 



1 



(49) 



and t = — (a + 1). In the region above the separatrix 
(X > Xa, corresponding to case B) the fluid behaves 
in a standard manner. The IS fixed point has repeller 
stability character and all trajectories are asymptotic to 
it in the past. The dS fixed point has attractor stability 
and is the generic future attractor. The KS fixed point 
has saddle stability in this region. In the region below 
the separatrix (X < Xa, corresponding to case C) the 
fluid behaves in a phantom manner. The KS fixed point 
has repeller stability character and the dS fixed point 
is an attractor. The M fixed point has saddle stability 
character. In general, phantom trajectories in this region 
are asymptotic to the Kasner singularity in the past and 
the generalized dc Sitter model in the future. 



B. Bianchi V cosmologies 

The models of type Bianchi V are a subclass of the 
Bianchi Class B models (a a ^ 0). These are homo- 
geneous and anisotropic cosmological models containing 



both the flat Friedmann (F) model and the open Fried- 
mann models (the Milne model for p = Q). In the case of 
Bianchi V models the group invariant orthonormal frame 
{u, e Q } can again be specialized further: 

ai = a 3 = n a p = 0, (50) 
oi + 0, (51) 
7 M a/3 = 0. (52) 

This in turn means that the spatial curvature scalar is: 



A R = -6af. 
Where the evolution equation for a\ is: 

di = —Ha\. 
The Einstein field equation (|17fl implies: 
<T a i = 0. 



(53) 



(54) 



(55) 



We will continue to use the normalized variables intro- 
duced in the previous subsection, x, y, z and r\. Addition- 
ally we introduce the new normalized curvature variable 
c: 

di 



The generalized Friedmann equation is now: 

v2 



y 



X 

3 



z 



(56) 



(57) 
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As with the previous case, this constraint equation can 
be used to reduce the dimension of the dynamical system: 
again we chose to remove y. The resulting set of evolution 
equation is then: 

x = — v/3x + 3z 2 + 9c 2 ((a + l)x + x 2 ) , (58) 



z 1 = -z^fZx + 3z 2 + 9c 2 , (59) 




(60) 



The fixed points of the new reduced system, the corre- 
sponding eigenvalues and their existence conditions (the 
energy density remains positive, x > 0, and all four vari- 
ables are real, x, y, z, c <E R) are given in Table ITTT1 The 
fixed points are identical to those found in the Bianchi I 
case. As with the Bianchi I case, in order to analyze the 
dynamics at infinity, we introduce the Poincare variables 
Xi. We then study the compactified and reduced system 
(A, Z, C), with the fixed points given in Table llVl 



TABLE III: Location, eigenvalues and existence conditions 
(x > and x, y, z, c g R) for the finite value fixed points 
of the reduced Bianchi V system. To simplify we use r = 
-(a + 1). 



The first two fixed points M and dS represent the 
Minkowski and generalized de Sitter points as in the 
Bianchi I case and retain the same stability character. 
The third point (CS) is the curvature singularity 3 , which 
is curvature dominated (C — > 1, or c — ► oo). The fourth 
point (IS) represents the isotropic singularity and the fi- 
nal fixed point (KS) represents the anisotropic Kasner 
singularity. We now discuss the state space diagrams for 
this system. The case a > — 1 is shown in Fig. OJa) (this 
corresponds to case A). The X — Z plane (C = 0) is 
the invariant manifold which represents the reduced and 
compactified Bianchi I state space, which is a manifold 
of the larger Bianchi V state space. The Z — C plane 
(X = 0) is the vacuum invariant manifold and the X — C 
plane (Z = 0) is the isotropic invariant manifold. The 
Minkowski point M is the future attractor for all trajec- 
tories in the region. The isotropic singularity (IS) is the 
generic past attractor, while the anisotropic Kasner sin- 
gularity (KS) and curvature singularity (CS) have saddle 
stability character. As in the Bianchi 1 case, in general 
the past singularity is isotropic. 

In the case a < — 1, depicted in Fig.Qtb), there are five 
fixed points, M, dS, CS, IS and KS. The state space is 
split into two regions by the separatrix manifold joining 
the points dS, CS and KS. This manifold is given by 
X = Aa, where: 



Name x, z, c 



Eigenvalues 



Existence 



M 



dS 



a G 



(0, 0, 0) 0, 0, 

(r, 0, 0) -tV3t, -V&r, -v/t-73 a < -1 



TABLE IV: Location and existence conditions (X > and 
X, Z, C G R) for the compactified and reduced Bianchi V 
system. To simplify we use r = —(a + 1). The * is shown 
when the point does not exist. 



Name 


X, z, c 


Existence Stability (a < — 1, a > — 


M 


(0, 0, 0) 


a e R 


Saddle, Attractor 


dS ^ 


75*- °< °) 


Q < -1 


Attractor, * 








CS 


(0, 0, 1) 


a e R 


Saddle, Saddle 


IS 


(1, 0, 0) 


a G R 


Repeller, Repeller 


KS 


(0, 1, 0) 


a G R 


Saddle, Saddle 



A A 



1 - C 2 - Z 2 



(61) 



and t = —(a + 1). In the Z — C plane (A = 0, vacuum 
invariant manifold) the dynamics are unaffected by the 
change in a. The A — C plane (Z = 0, isotropic invariant 
manifold) is divided into two regions by the separatrix. 
In the region above the separatrix (A > Aa, correspond- 
ing to case B) the fluid behaves in a standard manner. 
The IS fixed point has repeller stability character and 
all trajectories are asymptotic to it in the past. The dS 
fixed point is the generic future attractor. The CS and 
KS fixed points have saddle stability character in this 
region. In the region below the separatrix (A < Aa, 
corresponding to case C) the fluid behaves in a phan- 
tom manner. The KS fixed point is a repeller and the 
dS fixed point the attractor. The M and CS fixed points 
have saddle stability character. In general, trajectories 
are asymptotic to the Kasner singularity in the past and 
the generalized de Sitter model in the future. 



3 Strictly speaking, this points represents at the same time the 
true past curvature singularity of open Friedmann models (the 
left face of the sphere quarter in Fig.|3J as well as the coordinate 
singularity of the Milne model, the line X = Z = in the same 
figure. 
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FIG. 3: The compactified state space (X, Z, C) for the Bianchi V reduced system, including fixed points at infinity. These 
represent the isotropic singularity (IS), the Kasner anisotropic singularity (KS) and the curvature singularity (CS). (a) Left 
panel: The state space when a > — 1. (b) Right panel: the state space when a < — 1. It is clear that standard (non phantom) 
trajectories are asymptotic to the IS fixed point in the past. 



IV. DYNAMICS OF ANISOTROPY AND 
INHOMOGENEITY 

We have seen in the previous section that generic tra- 
jectories of Bianchi I and V anisotropic models with stan- 
dard (non-phantom) matter of type A and B asymptot- 
ically approach an isotropic singularity in the past, with 
the Hubble normalized shear vanishing, £ — > 0. Since 
the Bianchi I and V models are rather special and homo- 
geneous, we now want to see if this behavior is restricted 
to these models only, or if instead is stable against small 
generic inhomogeneuos perturbations. That is, we want 
to see if appropriately defined dimensionless quantities 
measuring the deviation from homogeneity in a covari- 
ant gauge invariant way decay to zero when the singu- 
larity is approached. Since we are only interested in the 
neighborhood of the singularity, we only need to consider 
perturbations of the spatially flat (K = 0) homogenous 
and isotropic cosmological models (the flat Friedmann 
models, i.e. the vertical X axis in Figs. and |3J) with 
quadratic EoS 



P 



ap 



£1 

Pc 



(62) 



Given this EoS, p + 3P > is always satisfied for high 
enough p, whatever a is, hence at high enough ener- 
gies the background evolution is non-inflationary, d < 0. 
Therefore we only need to consider large scale pertur- 
bations, i.e. scales much larger than the Hubble horizon 



-ff -1 , because for any given perturbation scale A there is 
always a time tn (the horizon crossing time) such that 
for t <C tn we have A » if -1 . 

For simplicity we now focus on models of type A only, 
with a > — 1; it can be shown that similar results can be 
obtained in case B, a < —1. For this case we re- write 
Eq. JHJ as 



p{a) 



where 



Pc(a + 1) 

(JL)3(a+l) _ l 



= j43(c«+1) 



(63) 



(64) 



is the finite value of the scale factor such that the energy 
density has a pole like singularity. A is a constant given 
by jnj and is expressed in terms of p , a , p c and a. 
Given that the background is a spatially flat Friedmann 
model the Friedmann equation (H 2 = p/3) gives directly 
the Hubble expansion scalar H, while a Q can be fixed 
arbitrarily, a freedom that we will use later. The sound 
speed for the fluid is given by: 



2 dP 2p 
c s = — = a+ — . 
op Pc 

The more commonly used EoS parameter w is: 



W 



P 
P 



Pc 



(65) 



(66) 
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The sound speed and EoS variable are then related by 
the following relationship: 



d = 2w 



(67) 



As said above, we will carry out our analysis in the long- 
wavelength limit A > H~ l , or k 2 /a 2 H 2 < 1 (A = a/k). 
This is analogous to dropping the Laplacian terms in the 
evolution and constraint equations. 



Thus in case A our parameter a could be replaced by 



the combination a = 2w - 
parameters w and cz. 



c s of the more commonly used 



B. Scalar perturbations 



A. Dimensionless variables and harmonics 

We use the covariant 1+3 approach j2^, H^, EIJ 
to study inhomogeneous and anisotropic gauge-invariant 
perturbations. In this approach, the information relat- 
ing to inhomogeneous and anisotropic perturbations is 
contained in a set of (1+3) covariant variables. These 
variables are defined with respect to a preferred time- 
like observer congruence u a . We use appropriate dimen- 
sionless density, expansion and curvature gradients which 
describe the scalar and vector parts of the corresponding 
perturbations: 



Z a 
C a 



-V a p, 
P 

Za - a V TT 



(68) 

(69) 
(70) 



The dimensionless expansion normalized quantities for 
the shear <7 ao , the vorticity uj a , the electric E ab and mag- 
netic H a b parts of the Weyl tensor are: 



H 



H 



a ab TI r W a E a i, 
— -TT J W a — TT 7 tab — -Jp 



"Hob = 



Hab 



where H is the Hubble scalar, H = a /a. 



H 2 ' 
(71) 

Here uj a = 



Vabc^ is t ne usual vorticity vector. In order to calculate 
all the relevant tensor perturbations it is useful to define a 
set of variables corresponding to the curls of the standard 
quantities. The curl of a given quantity will be denoted 
by an over-bar. 



S fl 6 = 77CurlS ab 
ti 



£ab = JjCUll£ ab 



Hab = -^cmYHab 



(72) 

The harmonics defined in [29J are used to expand the 
above tensors X a , in terms of scalar (S), vector (V) 
and tensor (T) harmonics Q. This results in a covariant 
and gauge invariant splitting of the evolution and con- 
straint equations for the above quantities into three sets 
of evolution and constraint equations for scalar, vector 
and tensor modes. The scalar, vector and tensor quanti- 
ties can be split such that: 



X 
X a 

X a b 



= X b Q 



Sr>s 



k- 2 X s Q s ab 



X v Q v a 



k- l X v Q v ab + X T Q T ab 



(73) 
(74) 
(75) 



The evolution and constraint equations are greatly 
simplified if we use the previously defined dimensionless 
time derivative: 



( )' 



l__d 

Hit 



da 



(76) 



The linearized evolution equations for scalar perturba- 
tions are: 



(77) 
(78) 
(79) 
(80) 



A' 


= 3wA- 


H ' 


Z' 


= -2Z - 


3H . 
~2~ 


E' 


= -3E- 


— (w+l)a, 


a' 


= -2(7- 


E 
H' 



The set of equations contain a mixture of Hubble normal- 
ized and non- normalized variables, however all solutions 
presented are appropriately normalized. In a addition we 
have a set of constraint equations for the scalar quanti- 
ties: 



C = -Aa 2 HZ + 2a 2 pA, 
3S = 2Z, 

e = a. 



(81) 
(82) 
(83) 



The system of equations is particularly difficult to solve in 
the above form. This is due to the form of p(a), H(a) and 
w(a) (the pole like behavior at a = a s ). In order to solve 
the system of equations we carry out a change of vari- 
ables. This is a parameter specific {p c and a) non-linear 
transformation of the scale factor. Once the solutions are 
obtained, they are again expressed in terms of the origi- 
nal normalized scale factor a. The solutions to the system 
of equations can be split into modes which are sourced 
by the gradient of the 3-curvature (C, which is constant 
in the long-wavelength limit) perturbation (subscript +) 
and those which are not (subscript — ). The solutions 
for the scalar contribution to the density and expansion 
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TABLE V: The large scale behavior of the scalar quantities. The constants of integration (Ai) correspond to the two 
independent modes. The non-zero constants Bi, Zi, Co, Si and £i can all be expressed in terms of Ai^, pc and a. The symbol 
5 represents the generalized hypergeometric function given in Eq. 1881 . The third column gives the behavior of the perturbation 
at the singularity. The last column gives the approximate late time behavior of the perturbation. The limits are given under 
the assumption that a > —1 and p c > 0. 

Quantity X- Mode X+ Mode Limit (a — ► a s ) Limit (a 3> o s ) 



A Ai Vpc(* + i)a3(. + i )-i A 2 M" +1 ^ ( " +1) - 1] (1 + ^5) 0,0 



V Pe(a + l)„3( n+ l)_i [p c (Q + l)q3(° + 1 )-l] , Bi E) 



£ £i ^ (a +y ti)-i g2 LPetn+ ija- 1] (1 + a.g) o,o 



a 2 






[p c (a + l)a 3 < Q 


+ 1) 


-i] 










Co 




[p c (a + l)a 3 ( Q 


+i) 


-i] 








[p c (a + l)<x 3 < a 
_i! 


+i) 


-i] 



z Zi V Pe (a+iy( «+l)-l ^ |,, ,. , ,.„ _ M(j_. i^£j () (J 

CO Co 0, Co 0, Co 



gradient's are of the form: 




The Aj's (i = 1,2) are the arbitrary constants of inte- 
gration and correspond to the two independent modes. 
The other constants Bi and Zi can be expressed in terms 
of the arbitrary constants Aj and the EoS parameters p c 
and a. The symbol S represents the generalized hyper- 
geometric function and is of the form: 



approximately: 

A_ » Aia 3(Q - 1)/2 , (90) 

A+ » A 2 a 3a+1 , (91) 

Z_ « Z ia 3(Q ~ 1)/2 , (92) 

Z+ m Z 2 a 3a+1 . (93) 

This is the corresponding behavior one would expect for a 
fluid with a linear EoS. The modes asymptotically evolve 
as there linear EoS counterparts at late times. The per- 
turbation modes for each of the scalar quantities can be 
found in Table [V] The constants of integration are the 
Aj's and all other constants {Zi, Ej, Si, Bi and Co ) can 
be expressed in terms of these constants and the EoS pa- 
rameters. The behavior of each mode is also given in the 
two limits, the first is at the singularity (a — + a s ) and the 
second is at late times (a 3> a s ). The limits are given 
assuming a > — 1 and p c > 0. 





'1 3a + 2 " 




"3" 


,1-b) 


( 


2' 3(a+l)_ 




2 



b = p c (a + l)a 3(Q+1) . (89) 

The generalized hypergeometric function is equal to unity 
at the singularity (2^1 ([i, j], [k], 0) = 1), where also 6=1 
after using the freedom in a to set A — [p c (l + a)] -1 . 
In addition the quantity (1 ± BjH/a) -* 1 as a > a s , 
the hypergeometric function is a correction to the solu- 
tions at the singularity and plays a subdominant role at 
late times. The perturbations all decay as we approach 
the singularity in the past (A± — > and Z± — > as 
a — > a s ). At late times (a 3> a s ) the various modes are 



C. Vector perturbations 

The linearized evolution equations for the vector com- 
ponents of the perturbations in the large-scale limit are: 



A' 


= 3wA 


H ^' 


(94) 


Z' 


= -2Z 


3H . 
2~ 


(95) 


oj' 


= (3c 2 - 


-2)w, 


(96) 


E' 


= -3E 




(97) 


a' 


= -2(7- 


E 

~ H' 


(98) 
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TABLE VI: The large scale behavior of the vector quantities. The constants of integration (Aj and Wi) correspond to the 
three independent modes. The non-zero constants Bi, Zi, Co, Hi, £j and £i can all be expressed in terms of Ai,2, Wi, p c and 
a. See the caption of Table IV! for further details. 



Quantity X- Mode X+ Mode Limit (a — > a 3 ) Limit (a ^> a 3 ) 



A 




A 2 


[p c (a + l)a 3 <' 


» + !)_!] 


(1 + 


2x3) 


0, 


8(a-l)/2 j ffl 3aH 


-1 






0- 




a } 




Z 


ry VPc(^ + l)a 3 ( a + 1 )-l 




[p c (a+l)o 3 (' 


<+D_l] 


(1- 




0, 


a 3(a-l)/2 j ffl 3aH 


-1 




7 — 




a 1—1 / 




c 









Co 






0, Co 


0, Co 




w 


117 lp c (c + l)a 3 ^+^-l) 5 / 2 

Wi aia + 8 















a (9a-l)/2 




H 


^ [p c ( a + l) a 3( = + l)-l] 















a 3(a-l) 




E 


Vpc(<i + 1)<i 3( ° + 1) -1 


E 2 


[p c (a+l)a 3 (< 


» + !)_!] 


(1- 


B 2 ^\ 


0, 


a 3(a-l)/2 j ffl 3cH 


-1 


ZjI 3 










£ 


c Vfc(a + l)a 3<a+1) -l 


£2 


[p c (a + l)a 3 < c 


. + !)_!] 


(1+ 




0, 


a 3(a-l)/ 2j fl 3cH 


-1 


61 3 






a / 





The constraint equations for the vector variables are 

C = XjrUZ + 2<r//A. (99) 



-Aa 2 HZ + 2a 2 pA, 
H = '^-{w + l)HW, 

3£ = 2Z, 
E = A. 



(100) 

(101) 
(102) 



The perturbation modes for each of the vector quantities 
can be found in Table IVll The constants of integration 
(A, and W\) correspond to the three independent modes. 
The non-zero constants Bi, Zi, Co, T~L\, Sj and Si can 
all be expressed in terms of the constants of integration 
and the EoS parameters. The last two columns give the 
behavior of the perturbation in the two limits. As in the 
case of the scalar contributions all perturbations decay at 
the singularity and asymptotically approach the standard 
linear EoS behavior at late times. 



D. Tensor perturbations 

The linearized evolution equations for the tensor com- 
ponents of the perturbations in the large-scale limit are: 



E = 


-3E 


a' = 


-2(7 


H' = 


-m 


E' = 


-3E 



3H . 
— {W 




(103) 


E 
H' 




(104) 


H 




(105) 


3H . 
— (W 


+ 1)H. 


(106) 



The constraint equations for the tensor variables are 

H = t, (107) 

U = 0. (108) 

The perturbation modes for each of the tensor quantities 
can be found in Table IVIII The constants of integration 
(Sj and Tii) correspond to the four independent modes. 
The non-zero constants Bi, £i, £i and Ej can all be ex- 
pressed in terms of the constants of integration and the 
EoS parameters. The perturbations all decay at the sin- 
gularity and evolve as there linear EoS counterparts at 
late times. 



V. CONCLUSIONS 

In this paper we have investigated the effects of 
a quadratic EoS in homogenous and inhomogeneous 
anisotropic cosmological models in general relativity, in 
order to understand if in this context the quadratic EoS 
can isotropize the universe at early times, when the initial 
singularity is approached. 

The first motivation for this work is that in the context 
of brane world cosmological models the term quadratic 
in the energy density that appears in the effective 4- 
dimensional equations of motion on the brane dominates 
at early times, driving the evolution to an initial isotropic 
singularity under reasonable assumptions. Therefore, un- 
der these assumptions, in the brane world scenario we 
do not have the classical isotropy problem: the isotropy 
we observe today is built in, that is it follows from 
generic initial conditions. This is at odds with the classi- 
cal and now widely accepted results in general relativity 
where general cosmological models with a lin- 
ear EoS are highly anisotropic in the past, approaching 
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TABLE VII: The large scale behavior of the tensor quantities. The constants of integration (Ei and Hi) correspond to the 
four independent modes. The non-zero constants Bi, Si, £i and £i can all be expressed in terms of £1,2, "Wi,2, pc and a. See 
the caption of Table IVI for further details. 

Quantity X- Mode X+ Mode Limit (a — > a s ) Limit (a ^> a s 

* a 13 a z V a / 



f £i Vp c (.+» 3 '^ l £2 L Pe ^-n;^ 11 ( 1 + £l H ) !U) 





+ l)a 3 < 


1 ( 


2 — 






. [Pc(a- 


fl)a 3( ' 


»+!)_!] ^ 








[Pc(<*H 


-l)a 3 < c 


.+!)_!]! 




a 3 




[Pc(«4 


l)a 3 <° 


+ D_1]2 




a» 






-l)a 3 C° 




ail 



H Hi [Ma+i)a 3(o+1) -i] ^ |J "'^ — ( I — — E ) ii. 1) 

g g - 1 fpc(^+i)^° +1 '-il g 2 M"+i)^° +1) -i] § (1 + ^5) n. i) 

g, Si lpc(a+1) ^ ( ° +1> ~ 11 v.. ■"■ ' 1 . ^-(1. - ~-5) U. 

H 



the singularity a-la BKL, i.e. with mixmaster dynamics, 
so that in general one needs special initial conditions to 
have isotropy at late times, even to initiate an inflation- 
ary phase [3 ■ 

Secondly, we also wanted to study the effects at early 
times and at high energies of the same type of EoS used 
in Paper I Q] , relaxingthe symmetries assumed there. 
The focus in Paper I [l| was on the possible use of a 
quadratic EoS as an effective way of representing a dark 
energy component or unified dark matter in the context 
of homogeneous isotropic models of relevance for the uni- 
verse at late times. Here we have used the simplified EoS 
P = ap + p 2 /pc, which still allows for an effective cos- 
mological constant and phantom behavior, and is general 
enough to analyze the dynamics at high energies. 

In summary, our aim here was to investigate if the be- 
havior of the anisotropy at the singularity found in the 
brane scenario is recreated in the general relativistic con- 
text if we consider a quadratic term in the EoS. That is 
to say, in the language of dynamical system theory, do 
we have a scenario where the generic asymptotic past at- 
tractor in general relativity is an isotropic model, thanks 
to a non- linear (specifically, quadratic) EoS. 

To this end, we have first summarized in Section [H] 
the classification of the energy density evolution in three 
classes, A, B and C, following from the energy conserva- 
tion for the EoS P = ap + p 2 / p c (a broader classification 
can be found in Section II of Paper I ID) . In Section ITTT1 
using dynamical system methods |26l l27j . we have stud- 
ied the Bianchi I and V cosmological models with a fluid 
with this EoS. We first considered the case when a > — 1 
(case A). The physically interesting quantity E, the Hub- 
ble normalized dimcnsionlcss measure of shear, describes 
the anisotropy. In the expanding Bianchi I and V mod- 
els this quantity always decreases at late times, as in the 
case with a linear EoS, approaching a Minkowski model. 



However, we can now have stages at early times in which 
E grows, vanishing at the initial singularity, and generic 
trajectories in the compactified state space asymptoti- 
cally approach in the past the isotropic fixed point IS, 
as depicted in Figs. [5] and |3] This behavior is similar 
to the quasi-isotropic expansion seen in the anisotropic 
brane-world models, but occurs for different reasons. The 
quasi-isotropic expansion is caused by the fact that the 
EoS becomes "super-stiff" at early times but decreases 
to finite values at late times. At the initial singular- 
ity (which is of Type III, see H2,|3^|) the anisotropy is 
always sub-dominant and we recover the standard linear 
EoS behavior at late times. In the case a < — 1, the state 
space can be divided in two regions. In the standard fluid 
region (X > X\, case B) the behavior at early times is 
similar to the a > — 1 case, with generic trajectories in 
the compactified state space approaching the isotropic 
fixed point IS of Figs. [21 and [3] At late times instead tra- 
jectories asymptotically approach generalized expanding 
de Sitter models. Finally, in the phantom fluid region 
(X < X\, case C) the anisotropy dominates at early 
times and the past singularity is anisotropic. Trajecto- 
ries are asymptotic to the expanding de Sitter model at 
late times. 

The submanifold S = in the state space for the 
Bianchi models represents homogeneous isotropic spa- 
tially flat Friedmann-Robertson- Walker models, with IS 
as past attractor. In order to see if IS can be regarded 
as the past attractor of a larger class of general models 
with no symmetries, we have studied in Section llVl using 
covariant and gauge invariant variables |28l . l29l y(J j the 
first order inhomogcneous anisotropic perturbations of 
these flat isotropic models, focusing on those of class A. 
The linearized equations have been solved in the large 
scale limit A ^> H^ 1 , neglecting the Laplacian opera- 
tors terms. As explained in Section^ only perturbations 
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larger then the Hubble horizon H 1 are relevant in our 
analysis. The main results of this section are given in 
Table's El EH and EH1 The tables contain the form of 
the solutions and asymptotic behavior of the physically 
relevant quantities for both early (a — > a s ) and late times 
(a 3> a s ). The quantities have been harmonically decom- 
posed into the standard scalar, vector and tensor com- 
ponents. These include the expansion normalized vortic- 
ity shear, and electric and magnetic contributions of the 
Weyl tensor. In addition, we consider the gradients of 
the energy density, expansion and three-curvature. As 
expected, at late times the perturbations evolve as if the 
fluid obeyed a linear EoS, i.e. the standard behavior is re- 
covered in this limit. The evolution of the perturbations 
at early times is our main result. Indeed, we find that all 
perturbations tend to zero near the singularity. Thus a 
perturbed quasi-isotropic model with generic initial con- 
dition evolves from the isotropic past attractor IS. This 
indicates that in general relativity with matter described 
by a non-linear EoS there are general cosmological mod- 
els, corresponding to a non zero measure subset of all 
possible generic initial conditions [l3 |. that isotropize as 
the initial singularity is approached, with IS as the past 
attractor. As for the brane models, this would be at odds 
with the standard BKL results [mlHHi]. 

Singularities are commonly regarded as signaling the 
failure of classical gravitational theories, requiring a 
quantum theory of gravity to deal with trans-Plankian 
energies. Even in effective theories, perhaps cosmo logy 
doesn't require them, as in the pre-big-bang scenario |34| . 
In any case, it is important to understand which is the 
typical classical dynamical behavior in the neighborhood 



of singularities, as this possibly signals what the generic 
outcome of a quantum theory should be at low (sub- 
Plankian) energies. In the BKL picture for fluids with 
linear EoS and other simple type of matter 0] this be- 
havior is mixmaster and the effects of matter become 
negligible close to the singularity: we may say that ef- 
fectively in this limit the space-time still tells matter 
how to move, but matter is no longer able to tell the 
space-time how to curve. Perhaps this behavior should 
be regarded as pathological, and peculiar of the energy- 
momentum tensor of linear EoS matter (or other sim- 
ple energy-momentum tensor) in general relativity. In 
practice we know close to nothing about matter fields at 
the highest (close to Plankian) cosmological energies, and 
these may require a non-linear effective EoS description 
in general, or maybe Einstein equations with an effective 
non-linear EoS are the Einstein frame version of some ef- 
fective low energy (sub-Plankian) theory. In both cases, 
our analysis of Einstein equations with a fluid with a 
quadratic EoS indicates that the typical initial state is 
isotropic. The singularity is matter dominated, and of 
a type that could admit an extension. Interestingly, the 
EoS considered here are a subclass of those used in Paper 
I 1J to model a dark energy component or unified dark 
matter, allowing for effective cosmological constants. 
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